Lifetime measurement of the metastable 3d 2 D 5 /2 state in the 
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We present a measurement of the lifetime of the metastable 3d 2 D 5 / 2 state in the 40 Ca + ion, using 
the so-called shelving technique on a string of five Doppler laser-cooled ions in a linear Paul trap. 
A detailed account of the data analysis is given, and systematic effects due to unwanted excitation 
processes and collisions with background gas atoms are discussed and estimated. From a total of 
6805 shelving events, we obtain a lifetime r = 1149 ± 14(stat.) ± 4(sys.) ms, a result which is in 
agreement with the most recent measurements. 
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I. INTRODUCTION 



The metastable 3d 2 D 5 / 2 state in 40 Ca + is interest- 
ing for such diverse areas of physics as atomic structure 
calculations, optical frequency standards, quantum infor- 
mation and astronomy. For atomic structure calculations 
it is an important test case for the study of valence- 
core interactions and core-polarization effects 1]. The 
long lifetime of the 2 D 5 / 2 state implies a sub- Hz nat- 
ural linewidth of the 729 nm electric quadrupole transi- 
tion to the 2 Si/2 ground state, thus making this tran- 
sition an attractive candidate for an optical frequency 
standard 0. 01 In addition, its long lifetime makes the 
2 D$/2 state a suitable choice for one of the two states 
in a quantum bit in connection with quantum computa- 
tion Q. Finally, in astronomy, the 2 D 5 /^ state has been 
used in a study of the (3 pictoris disk 5], and it is also 
commonly used in the study of Seyfert 1 galaxies and T 
Tauri stars 0. As a consequence, the natural lifetime of 
the 2 D 5 /2 state has attracted much attention in recent 
years, but unfortunately the results of the many mea- 
surements [lJ 3 jpj llHinOOlF4]| and theoretical 
calculations 0, |g, [Tj, [T|| E3, E 03 of the lifetime are 
scattered over a rather broad range (see Fig. 1 in Ref. 
for an overview). Among the experimental results, the 
lifetime found by Barton et al. 7], using the so-called 
shelving technique on a single trapped and laser-cooled 
i0 Ca + ion, has the smallest error bars, with an estimated 
lifetime of r = 1168 ± 7 ms. 

In this paper we present a measurement of the lifetime 
of the 3d 2 D 5 / 2 state, using the same shelving technique 
but on a string of five ions. The measurement results in 
a lifetime of r = 1149 ± 14(stat.) ± 4(sys.) ms. 

The experimental setup and the experimental proce- 
dure are described in Section [H] In Section ITTT1 we give 
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a detailed account of the data analysis, including a de- 
scription of the maximum likelihood method used for the 
statistical data analysis and a discussion of systematic 
errors with emphasis on radiation effects and collision ef- 
fects. Finally, in Scction llVl the estimated lifetime based 
on the measurements is given and discussed. 



II. EXPERIMENTAL SETUP 

A sketch of our experimental setup is shown in Fig. ^ 
The linear Paul trap used in the experiments is situated 
in a stainless steel ultra-high- vacuum chamber. The trap 
consists of four cylindrical gold-coated stainless steel rods 
arranged in a quadrupole configuration, where two di- 
agonally opposite rods are separated by 7.00 mm. The 
rods are 8.00 mm in diameter, and each rod is sectioned 
into three parts, where the center piece is 5.40 mm long, 
and the two end pieces are 20.00 mm long. By applying 
an RF-voltage to two diagonally opposite electrode rods 
and applying the same voltage with opposite phase to 
the other two electrode rods, we obtain an effectively 
harmonic radially confining potential. In the experi- 
ments reported here, the peak-peak amplitude of the RF- 
voltage is 600 V and the frequency is 3.894 MHz, which 
results in a radial trap frequency of ui r w 2-7T • 550 kHz. 
Axial confinement is provided by a 580 mV DC- voltage 
on all eight end electrodes, yielding an axial trap fre- 
quency u> z w 2-7T • 50 kHz. By adding additional small 
DC-voltages to the individual electrodes, we can care- 
fully center ions in the trap. At the center of the trap 
a collimated atomic Ca beam is crossed by a photoion- 
izing laser beam, and 40 Ca + ions are produced by an 
isotope-selective resonance-enhanced two-photon ioniza- 
tion process [20l l2lj and loaded into the trap. The Ca 
atomic beam originates from a calcium sample contained 
in an oven. The Ca atoms effuse out of a hole in the oven, 
and by a number of skimmers a collimated atomic beam 
is produced. An oven shutter enables us to fully block 
the atomic beam. In the present experiments the oven 
was heated to 420° C, under which conditions the cham- 
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ber pressure was measured by an ion gauge 27] to be 
6.0 • 10~ n Torr. After having loaded five ions and forced 
them onto a string by adjusting the trap parameters and 
applying Doppler laser-cooling as described below, the 
oven temperature was reduced, and the oven shutter was 
closed. During a one hour measuring session, the pres- 
sure dropped to about 3.6 • 10~ n Torr. The ions are 
Doppler laser-cooled on the 4s 2 Si/2-4p 2 -Pi/2 transition 
with typically 15 mW of 397 nm light from a frequency- 
doubled TkSapphire laser. From the Ap 2 P 1 / 2 state, the 
ions will not always decay back to the 4s 2 S]/ 2 state 
but sometimes decay by a dipole-allowed transition to 
the 3d 2 D 3 / 2 state, from where they are repumped back 
to the ip 2 Pi/ 2 state using a diode laser at 866nm, see 
Fig. [21 The ions are observed by collecting fluorescence 
light, emitted at 397 nm during the Doppler cooling pro- 
cess, with a Nikon objective lens for MM40/60 measuring 
microscopes (lOx magnification, f-number^ 1.7) placed 
about 5 cm above the trap center (outside the vacuum 
chamber). The collected light is amplified by an image 
intensifier and imaged onto a CCD-camera, resulting in 
an all over magnification of 13.5 for the entire imaging 
system |28|. The CCD-camera has a 50 ms exposure time, 
and digitial images of the ions with 12-bit resolution are 
recorded to the RAM of a personal computer at a frame- 
rate of 17.610 Hz. In order to measure the lifetime of 
the 3d 2 Z3 5 /2 state, during the Doppler cooling process 
we excite ions in the 3d 2 D 3 / 2 state to the Ap 2 P 3 / 2 state 
using a diode laser at 850 nm, as shown in Fig. |2] From 
the 4p 2 P 3 / 2 state the ions can spontaneously decay to 
the 3d 2 D 5 / 2 state. When the ion is in the 3cZ 2 -D 5 / 2 
state, the optically active electron of the ion is said to 
be shelved, and the fluorescence on the 397 nm transition 
is quenched. The time spent by the ion in the 3d 2 D 5 / 2 
state before it decays back to the ground state, will in 
the following be called a shelving period. By continuously 
applying the two cooling lasers and the shelving laser at 
850 nm, we obtain characteristic fluorescence signals like 
the one shown in Fig. [3] 
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FIG. 1: Experimental setup, see text for details. M: mirror, 
DM: dichroic mirror, PBS: polarizing beamsplitter, A/2: half- 
wave plate. 




FIG. 2: Relevant levels and transitions in 40 Ca + . Doppler 
laser cooling is performed with lasers at 397 nm and 866 nm. 
A laser at 850 nm excites the ion from the 3d 2 D3/2 state to the 
4p 2 P 3 / 2 state, from where it can decay to the 3d 2 D 5 / 2 state. 
The dashed line indicates the electric quadrupole transition 
by which the ion can decay from the 2 D 5 / 2 state to the ground 
state. 



The lifetime of the 2 D§/ 2 state can easily be estimated 
from the distribution of the shelving periods, since we 
expect it to be exponential with a time constant equal to 
the lifetime t. 

The duration of a shelving period can be determined 
by dividing the number of consecutive frames, where an 
ion is shelved, by the frame-rate of the CCD-camera. 
Our raw data are digital images of five ions, as shown 
in Fig. [31a, and our main dataset consists of ~ 200.000 
such images taken in three one hour experimental runs. 
The details of obtaining the distribution of the shelving 
periods from these images is described below. 



From the images in Fig. [21 a, it is evident that the ions 
are spatially well resolved, and that a 'region of interest' 
(ROI) around each ion can be defined. The ROI is 13 x 13 
pixels, corresponding to a region of 9.5 /imx9.5 /im in the 
trap region. In order to establish a fluorescence data- 
point which reflects the real ion fluorescence rate within 
a given image frame, we simply integrate the pixel values 
within the ROI. This is a valid measure, since the image 
intensifier and the CCD-chip have a linear response to 
the fluorescence collected by the objective lens. In the 
following the ion fluorescence will be given in integrated 
pixel values (IPV). 
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FIG. 3: (a) A sequence of CCD-images including a shelv- 
ing event for the central ion. The white square indicates the 
'region of interest' (ROI) within which the fluorescence is in- 
tegrated to obtain the signal shown in (b). (b) Fluorescence 
signal from the central ion for an image sequence containing 
the nine images above. The datapoints corresponding to the 
nine images are indicated by open circles. The dashed line 
indicates the threshold level discussed in the text. 



As shown in Fig. the distribution of the fluorescence 
datapoints from a single ion is characterized by a rather 
sharp peak at a low fluorescence level and a Gaussian 
distribution of datapoints around a level of ~ 90000 IPV, 
originating from cases where the ion is scattering 397 nm 
light during a whole frame. In the following we define 
the fluorescing level as the center of the Gaussian distri- 
bution. The standard deviation of the Gaussian distri- 
bution is a ~ 7500 IPV, which is mainly set by image- 
intensifier noise, but also by laser intensity and frequency 
drifts during the experimental run, and by the finite ion 
temperature. The peak at the low fluorescence level, or 
the background level, is due to datapoints originating 
from cases where the ion is shelved during a whole frame. 
This background peak is asymmetric with a sharp left 
edge, since there is essentially no 397 nm light scattered 
into the ROTs, when the ions are shelved. 




20000 40000 60000 80000 100000 120000 

Integrated fluorescence [IPV] 

FIG. 4: Distribution of datapoints from a single ion, binned 
in intervals of 500 IPV. The maximum of the narrow back- 
ground peak near 4000 IPV is not shown but has a value of 
8000. From the Gaussian distribution around 90000 IPV with 
a width of a ~ 7500 IPV, we can define a fluorescing level as 
the center of this distribution. The inset shows the interme- 
diate points between 0IPV and 80000 IPV. 



In order to establish a distribution of the shelving pe- 
riods, we need to introduce a threshold level to discrim- 
inate between fluorescence datapoints corresponding to 
frames where the ion is fluorescing (above the thresh- 
old level) or shelved (below the threshold level). From 
the distribution of fluorescence datapoints in Fig. 0] it 
is evident that we can choose a threshold level which 
can be used to unambiguously distinguish between the 
background level and the fluorescing level of an ion. The 
specific choice of threshold level is discussed in the fol- 
lowing. 

From the inset of Fig. 21 it can be seen that the dis- 
tribution of fluorescence datapoints contains some points 
with a value of the integrated fluorescence between the 
background level and the fluorescing level. These inter- 
mediate datapoints arise from images where the ion is 
only fluorescing during a fraction of the exposure time 
of the CCD-chip. This occurs naturally when an ion is 
shelved or decays to the ground state during the exposure 
of a single frame. Since a particular choice of thresh- 
old level decides whether an intermediate datapoint is 
counted as belonging to a shelving period or not, it af- 
fects the precise distribution of the shelving periods. For- 
tunately, the choice of another threshold level on average 
only adds a constant amount to all the shelving periods, 
and therefore the decay rate extracted from their distri- 
bution is not influenced by the choice of threshold level. 
This fact allows us to choose a threshold level in a broad 
range between the background level and the fluorescing 
level. 

Unfortunately, intermediate fluorescence datapoints 
also occur in the following three situations: i) shelv- 
ing followed by fast decay to the ground state, ii) decay 
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to the ground state followed by fast re-shelving, iii) a 
shelved and an unshelved ion change places, e.g., due to 
a weak collision with a background gas atom or molecule 
or due to the finite temperature of the ions. Examples of 
these events, all taking place within one or two camera 
frames, are shown in Fig0a-c. Hence, before we choose 
the threshold level to be used for extracting the distribu- 
tion of shelving periods, we have to consider these three 
types of events in some detail. 

i) When an ion is shelved at a certain instant of a frame 
and decays back to the ground state within that same 
frame or the next, then the signal does not necessarily 
fall below the threshold level, as the example in Fig.[5]a. 
shows. Hence these events may not even be counted as 
shelving events. In the data analysis we simply account 
for such events by discarding all periods with a duration 
of one or two frames from the shelving period distribu- 
tion. Since the distribution is expected to be exponential, 
we can subsequently displace the remaining distribution, 
so that periods originally having a duration of n frames 
(n > 3) are set to be periods with a duration of n — 2 
frames. 

ii) If an ion decays and is quickly re-shelved, we can- 
not be confident that the signal rises above the threshold 
level, see Fig. EJb, and thus two shelving periods can 
appear as one longer shelving period, which would arti- 
ficially increase the extracted lifetime. Noting that such 
events require a high shelving rate, it is possible to ex- 
clude them on probability grounds. First, we set a low 
threshold level, T = 0.1 • (F~B) + B, where F is the fluo- 
rescing level, and B is the background level, which means 
that only the very fastest re-shelving events do not rise 
above the threshold level. Second, we require that after a 
decay, which ends a shelving period, there must be twenty 
consecutive datapoints above the threshold level; other- 
wise it is not counted as a shelving period in the data 
analysis. This requirement is only likely to be fulfilled 
with a low shelving rate, thus reducing the risk of ac- 
cepting shelving periods where a quick re-shelving event 
has happened. The probability for an event of decay and 
quick reshelving, with the signal not rising above the se- 
lected threshold level, followed by a decay and twenty 
consecutive datapoints above the threshold level, is be- 
low 2 permille regardless of the shelving rate. Hence this 
method gives at most a systematic error of -2 ms to the 
final result, and we will take the systematic error to be 
— 1 ± 1 ms. 

iii) As mentioned, two ions may change place, e.g., 
due to a weak collision with a background gas atom or 
molecule. Fig. 0c shows how a shelved and an unshelved 
ion changing place, effectively cut one long shelving pe- 
riod into two shorter ones, which artificially shortens the 
extracted lifetime. Therefore, in such events we restore 
the position of the ions to obtain a single shelving event, 
as shown in Fig. |5Jd. Positions are only restored if the 
value of the fluorescence datapoints for the two ions in- 
volved adds up to the fluorescing level F within ±2a of 
this level. We find about 200 such events in our main 



dataset. Since statistically the fluorescence from the two 
ions does not add up to the fluorescing level within ±2<r 
in all events, we estimate a systematic error to the life- 
time of +2 ms. True events of one ion shelving and an- 
other decaying within the same fraction of a frame will, 
however, occur, and erroneously be corrected by this pro- 
cedure. We estimate a systematic error of —10 ms due to 
the erroneously corrected events. All together the sys- 
tematic error due to ions switching place is —8 ms, with 
an estimated uncertainty of ±4 ms. 

During data-taking it happened that the ions heated 
up, so the ion-string became unstable. Such events are 
clearly visible on the images of the ions and were cut out 
of the dataset before performing the data reduction pro- 
cess described above. Likewise, periods where the lasers 
were adjusted or unstable are not considered in the fur- 
ther data analysis. 

After the data reduction process described above and 
using a threshold level of T = 0.1- (F— B)+B, we extract 
a distribution of the shelving periods for each ion, binned 
in time intervals of At = 56.786 ms, which is the inverse 
of the frame-rate. The counts in equivalent bins for all 
the ions in the three experimental runs are then added to 
obtain a single distribution, from which a decay rate can 
be determined. This distribution is shown in a histogram 
in Fig. HJ with the bins shifted such that the histogram 
has its origin at time zero. 



B. Statistical analysis 

From the histogram in Fig. we can infer the lifetime 
r of the 2 D 5 /2 state by assuming an exponential distribu- 
tion with a decay rate given by the inverse lifetime of the 
2 Z? 5 /2 state. The histogram is comprised of 6805 shelving 
events with duration up to more than 11 seconds. Since, 
at large times, there are only a small number of events in 
the distribution, a least squares fitting method is inap- 
propriate, and hence we employ a maximum likelihood 
estimate instead. 

The probability, pi, of falling into column i in the his- 
togram, i.e., having a shelving period of duration t, which 
fulfills ti < t < ii+i, where ti = i ■ At, is 

Pi = l -e- l l T dt = e- u ' T (l - e- A */ r ) . (1) 

Mathematically, the p^s define a proper probability dis- 
tribution, since ^2 i=0 Pi = 1- The sum starts at i = 0, 
corresponding to the origin of the distribution in the his- 
togram. Since the sum extends to infinity, we must, in 
principle, be able to measure infinitely long shelving pe- 
riods. In the case of real experiments, where one is only 
able to measure periods up to a finite duration, T max , 
the pi's should formally be renormalized by the factor 
[1 — exp(— Tmoj/r)] -1 . In our case where T max ~ 1 hour 
and r ~ 1 s we can, however, safely neglect this factor. 
In order to make the maximum likelihood estimate, we 
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FIG. 5: Examples of events which are problematic in connection with the data analysis. In all graphs the dashed line indicates 
the chosen threshold level of ~ 12500 IPV. (a) Shelving and fast decay; corresponding frames indicated by open squares, (b) 
Decay and fast re-shelving; corresponding frames indicated by open squares. Here the signal does rise above the threshold level, 
but for an even faster event this may not be the case, (c) A shelved and an unshelved ion change places, (d) The signals in (c) 
with the positions of the ions restored, i.e., the signals from Ion 1 and Ion 2 are interchanged after the crossing. At the frame 
where the two ions change places and the preceding frame, the fluorescence signal from Ion 1 and Ion 2 are given artificial 
values of 120000 IPV and IPV. 



now introduce the likelihood function 

L = m y\p± 



.i 



(2) 



i=0 



where rii is the number of shelving events in the i'th 
column of the histogram, and N = X^o n i IS ^ ne total 
number of shelving events. By maximizing L (or rather 
ln_L) with respect to r, we find 



At 



i)' 



(3) 



where < t > — ^njtj/iV is the mean duration of the 
shelving periods. The variance of the lifetime is 

-l 

/ M In T. \ 

Var(-r) 




Using Eq. © and Eq. J3J, the lifetime and the statisti- 
cal uncertainty can be determined solely from At, < t > 
and N, and we find t = 1154 ± 14 ms. The exponential 
distribution based on the maximum likelihood estimate 
is plotted as a solid line on top of the histogram in Fig. El 

As a check of the efficiency and validity of our data 
reduction process, we have performed a test of goodness- 
of-fit. To find the goodness-of-fit we employ the Kol- 
mogorov test, which is based on the Empirical Distribu- 
tion Function. A gen eral overview of EDF statistics can 
be found in Ref. [23, which also gives tables of signif- 
icance levels appropriate for cases like ours, where the 
lifetime is estimated from the data set. We find from our 
data a value for the Kolmogorov test of 0.44, clearly be- 
low the 10% significance level value of 0.995, thus showing 
that the fit is good. We note for completeness that our 
data set is binned, whereas the values in Ref. [2^| are for 
continuous data. However, since we have more than 100 
channels this is not expected to play any significant role. 
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FIG. 6: Histogram over the 6805 shelving events obtained 
after data reduction. Note that the number of shelving events 
is on a logarithmic scale. The error bars are the squareroot of 
the number of shelving events. The solid line is the maximum 
likelihood estimate of r = 1154 ms. 



C. Radiation, collisions and other potential 
systematic effects 

Apart from the systematic errors concerning the data 
reduction, already discussed in Section lill Al there are a 
few other relevant systematic errors, which will be dis- 
cussed in this section. 



Radiation 



From the level scheme of Fig. |21 it is obvious that any 
radiation which couples the 2 D 5 / 2 state to the 2 P3/2 
state may deplete the 2 D 5 / 2 state by excitation to the 
2 P 3 / 2 state followed by decay to the 2 5 , 1 / 2 ground state 



or the 2 D s /2 state and result in a measured lifetime 
shorter than the natural lifetime. As a consequence, the 
occurrence of such radiation must be considered and, if 
possible, reduced. Blackbody radiation at the relevant 
transition wavelength is negligible at room temperature, 
as also discussed by Barton et al. 0. Other 'thermal' 
sources such as the ion gauge, roomlight and computer 
screens were off during the measuring sessions, except for 
the screen where we watched the images of the ions. This 
screen was facing away from the vacuum chamber and is 
hence not expected to cause any problems. Another class 
of influencing light sources is the lasers used in the experi- 
ment. In particular the broad background of the emission 
spectrum of the 866 nm and the 850 nm diode lasers con- 
tains 854 nm light resonant with the 2 D 5 / 2 ~ 2 Pz/2 transi- 
tion. This radiation source was first recognized by Block 
et al. , and as noted by Barton et al. the discrepancy 
between their own and many of the earlier measurements 



is possibly due to this previously unrecognized source of 
error. To reduce the level of 854 nm light from the diode 
lasers as much as possible, a long-pass filter was inserted 
in the 866 nm bcamlinc and adjusted (by changing the 
angle of incidence) to 30% transmission at 866nm and 
< 5 • 10~ 4 at 854nm, and a grating (1200 lines/mm) 
with a spectrally selective diaphragm (0.9 mm diameter, 
1080 mm from the grating) was inserted in the 850 nm 
beamline. The filter, the grating and the diaphragm are 
shown in Fig. ^ The overall power of the lasers was re- 
duced to ~ 18 nW and ~ 1.7 mW at the place of the ions, 
for the shelving laser and the repumping laser, respec- 
tively. The waist size of the beams was 670 /im x 700 /im 
for the repumping laser and 360 fim x 430 /im for the 
shelving laser. 

In order to estimate the lifetime reduction due to ra- 
diation emitted from the diode lasers, we calculate the 
rate of de-shelving from the 2 D 5 / 2 state, i.e., the excita- 
tion rate from the 2 £>5/2 state to the 2 i"3/2 state times 
(1 — b), where b = 0.068 is the branching ratio for decay 
from the 2 i"3/2 state back to the 2 D 5 / 2 state [24). 

First we consider de-shelving due to the repumping 
laser, which we split into two contributions, one from 
the off-resonant 866 nm radiation and one from the near- 
resonant background radiation around 854 nm. The de- 
shelving rate due to the 866 nm radiation was calculated 
by Barton et. al. to J 866 • 9.4 • 10" 5 s^/froWnim" 2 ), 
where Isee is the intensity of 866 nm radiation. With 
Is66 = 2.3mWmm~ 2 we find a de-shelving rate of 2.2 • 
10 _4 s _1 , yielding a lifetime reduction of 0.3 ms. 

In order to calculate the excitation rate from the 2 -D 5 / 2 
state to the 2 P3/2 state due to radiation near 854 nm, 
we first consider the rate at a given frequency ojl, as 
expressed by Eq. (2) in Ref. 0: 



2J 2 + 1 7T 2 C 3 . I . 

Rl2 = 2J^ih^- 2 Ml -c 9{uJL 



(5) 



where Ji = 5/2 and J 2 = 3/2 are the total angular mo- 
menta of the involved levels, uji 2 = 27rc/854.209 nm [24| is 
the transition frequency, A 2 x = 7.7 • 10 6 s _1 the Ein- 
stein coefficient for spontaneous decay from the 2 i"3/2 
state to the 2 D 5 / 2 state, I the intensity of the incom- 
ing radiation, and g(u>L — ^12) a normalized Lorentz- 
distribution. Assuming a flat background spectrum of 
the diodelasers, i.e., the intensity per frequency interval 
AI/Au>l is constant, we can integrate Eq. JSJ) over fre- 
quency u>l and find an excitation rate of 



R = 



2J 2 



2 Ji + 1 hu 



1 7T 2 C 2 , Al 

-A 21 



12 



Auj l 



(6) 



Using a diffraction grating we have measured 
AI/Au> L < 0.19nW/(mm 2 • GHz) near 854 nm at the 
power used in the experiment, and we then find R — 
0.77 s -1 . Multiplying this rate by 1 — b and the trans- 
mission of 5 ■ 10~ 4 of the long-pass filter, we find a de- 



shelving rate of 3.6- 10 
of 0.5 ms. 



! s 1 , yielding a lifetime reduction 
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As a check of this estimate, we performed an experi- 
ment, again with five ions, with the cooling lasers on but 
without the shelving laser. Even without the shelving 
laser the ions may be shelved due to radiation from the 
866 nm rcpumping laser, which couples the 2 D 3 / 2 state 
to the 2 P3/2 state. From the observed shelving rate, we 
can then find the excitation rate on the 850 nm 2 D 3 / 2 - 
2 P 3 /2 transition and compare it to the calculated excita- 
tion rate for the 854 nm 2 D 5 / 2 - 2 P3/2 transition, taking 
the different linestrengths into account. In this exper- 
iment the power of the 866 nm laser was 7.3 mW, the 
waist was as above, the long-pass filter was removed and 
the oven-shutter was open. In about 35 minutes we ob- 
served 12 shelving events, i.e., the observed shelving-rate 
is 5.8 ■ 10~ 3 s -1 . Taking into account the number of ions, 
the population of the D-state (~ 1/3, since both cool- 
ing transitions are saturated) and the branching ratio &, 
we find the excitation rate on the 2 D 3 / 2 - 2 P 3 /2 transi- 
tion: (3/56)- 5.8 •10" 3 s- 1 =5.1-l(r 2 s- 1 . By multiply- 
ing this number with the relative linestrength between 
the 2 D 5 / 2 - 2 P 3 / 2 transition and the 2 D 3 / 2 - 2 P 3 / 2 transi- 
tion of 8.92 |2J|, we find an estimated value of 0.45 s _1 



1,00- 



for the excitation rate on the D, 



P 3 /2 transition, 



5/2" -f3/2 

which should be compared to R = 0.77 s _1 calculated 
above. The two numbers are not expected to be equal 
but only of the same order of magnitude, since the tran- 
sition wavelengths are different, and the diode laser back- 
ground is not necessarily equally strong near 850 nm and 
854 nm. Also some of the shelving events may be due to 
collisions, as discussed below. Nevertheless, the numbers 
agree within a factor of 2 and our calculated estimate of 
the excitation rate, and hence the de-shelving rate seems 
to be reasonable. 

Above we considered de-shelving due to the repump- 
ing laser. In exactly the same way we could consider 
de-shelving due to the shelving laser. However, since 
the intensity of the shelving laser is much lower than 
for the repumping laser, and the grating and the di- 
aphragm strongly reduce the level of 854 nm light, the ex- 
citation rate is expected to be extremely small. To check 
this we performed additional lifetime measurements at 
three different power levels of the shelving laser, 18 nW, 
193 nW and 2081 nW, yielding lifetimes of 1146 ± 24 ms, 
1160 ± 29 ms and 1092 ± 27 ms, respectively. When using 
higher power, the laser was detuned from resonance to 
get a similar shelving rate in all experiments. In Fig. 
the linear fit to the decay rates corresponding to the life- 
times shows that there is a weak dependence of the de- 
cay rate (or lifetime) on the 850 nm power. However, 
with only 18 nW the lifetime is only reduced by 0.5 ms. 
The results of the two low power measurements and the 
crossing at zero power are in agreement with the result 
obtained from our main dataset. The measurement at 
18 nW was in fact performed at the same shelving laser 
power as the measurements for the main dataset. 

All together, we include a systematic error of +1± 1 ms 
in our final result due to de-shelving from the two diode 



O 
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FIG. 7: Decay rate measurements at three different power 
levels of the shelving laser. A weighted least-squares linear 
fit to the measured decay rates yields 1/r = 0.866(8) s _1 + 
0.023(7) s-VnW-Psso. 



lasers. 



2. Collisions 

Another systematic effect which shortens the mea- 
sured lifetime is collisions with background gas atoms 
and molecules. There are two relevant types of colli- 
sions: fine-structure changing (j-mixing) collisions and 
quenching collisions. In a j-mixing collision the inter- 

3/2 



nal state can change from the D 5 / 2 state to the D 
state, or vice versa. In a quenching collision the internal 
state changes from the 2 D 5 / 2 state to the 2 6*1/2 ground 
state. Both types of collisions deplete the 2 D 5 / 2 state 
and hence shorten the measured lifetime, j-mixing and 
quenching rate constants (Tj and To) in the presence 
of different gases are given in Ref. [25| and references 
therein. Quite generally j-mixing collisions are found to 
be an order of magnitude stronger than quenching col- 
lisions. From a restgas analysis 29], we know that the 
restgas in our vacuum chamber is mainly composed of 
H2 and gas of 28 atomic mass units, i.e., N2 or CO. In 
the restgas analysis it was not possible to distinguish 
between N2 and CO, since the chamber pressure was 
so low that the signal from the atomic constituents of 
these molecules could not be observed. Knoop et al. [2{| 
found the following rate constants in units of cm s for 

r 3 -(H 2 



collisions with H 2 and N 2 : r 3 -(H 2 ) = (3 ± 2.2) • 10~ 10 , 
r Q (H 2 ) = (37 ±14) • 10- 12 , Fj(N 2 ) = (12.6 ± 10) ■ 10~ 10 
andr Q (N 2 ) = (170±20) • 10" 12 . The measurements were 
performed on a cloud of relatively hot ions, as compared 
to the laser cooled ions considered in this paper. As 
noted in Ref. pEj , other measurements of j-mixing with 
H2 at different collision energies give similar results, so 
we may expect only a weak energy- dependence for the j- 
mixing collisions, and therefore we use the values given in 
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Ref. |25j . On the other hand, higher quenching rates are 
expected at lower temperatures [2^, but to our knowl- 
edge there are no measurements of that for Ca + . Unfor- 
tunately, we are not aware of any similar measurements 
with CO, and in our estimate of the collision rate below, 
we therefore assume that the mass 28 restgas is N2. The 
relatively large rate constants of N2 found in Ref. [25| 
indicate that at least this assumption probably does not 
lead to a large underestimate of the collision induced de- 
shelving rate. 

Assuming that the values given in Ref. |25l | are applica- 
ble, we can estimate the collision-induced lifetime reduc- 
tion. At a pressure of 5 • 10~ n Torr, taking into account 
the sensitivity to different gases of the ion gauge and us- 
ing the result of the restgas analysis, we find that the 
restgas is composed of 54%H2 and 46%N2, which yields 
a total collision induced de-shelving rate of 2.3 • 1CP 3 s^ 1 , 
and a systematic error to the lifetime of +3 ms with an 
estimated uncertainty of ± 1 ms . 

From the measurement of 12 shelving events in 35 
minutes without the shelving laser on, as described in 
Sec. IIII C II above, we can obtain an upper limit for the 
j-mixing collision rate if we assume that all the observed 
shelving events are due to j-mixing collisions inducing 
a transition from the 2 D 3 / 2 state to the 2 D 5 / 2 state 
(with rate 735). Again taking the number of ions and 
the population of the 2 D 3 / 2 state into account, the col- 
lision induced shelving rate is 3.5 • 10 _3 s _1 . Since the 
oven shutter was open in that experiment, thus allowing 
collisions with neutral calcium atoms as well, we expect 
this number to be larger than under the conditions for 
the 'real' lifetime measurements. The transition rate, 
753, for the 2 I?5/2- 2 £>3/2 de-shelving transition is ex- 
pected to be given by 753 = 2735/3, owing to the principle 
of detailed balance, so the collision induced de-shelving 
rate is 2.3 • 10 s . This is our upper limit for the j- 
mixing collision rate, which should be compared to our 
estimate above of the total collision-induced de-shelving 
rate, which is dominated by contributions from j-mixing 
collisions. Somewhat fortuitously, the numbers are equal, 
and hence our calculated estimate seems to be reason- 
able. 



3. Other effects 

When observing shelving events from a string of ions, 
one might consider if there are any correlations in the 
decay of the individual ions from the 2 D 5 / 2 state to the 
ground state, which could influence the measured life- 
time. In the experiments by Block et al. 8] indications 
of correlated decays from the 2 D 5 / 2 state were observed, 
manifested as an overrepresentation of events where sev- 
eral ions decay at the same time. A later detailed ex- 
periment by Donald et al. fl4| showed, however, no such 
correlations. Apart from sudden bursts of 854 nm ra- 
diation, the only reasonable physical mechanism which 
could lead to correlations is so-called subradiant and su- 



perradiant spontaneous emission due to interference in 
the spontaneous decay of two or more ions |2fij . In 
the simple case of two ions, superradiant and subradi- 
ant spontaneous emission is characterized by a relative 
change in the normal (single ion) decay rate of the order 
of ±sm(kR)/kR (when kR > 10), where R is the ion-ion 
distance, k = 2n/X with A = 729 nm in our case j^D]. For 
more ions the effect is of the same order of magnitude. In 
our case 1/kR w 6 • 1CP 3 so just from this argument the 
effect is small, but not negligible. The interference effect, 
however, relies on creating and maintaining a superposi- 
tion state of the form |±) = (\S 1 D 2 )±\D 1 S 2 ))/V2, where 
S and D indicate the internal state, 2 Si/ 2 or 2 D 5 / 2l and 
indices 1 and 2 relate to Ion 1 and Ion 2. In our ex- 
periment such a superposition state can only be created 
by a random process since the D 5 / 2 state is only pop- 
ulated through spontaneous emission, and consequently 
the relative phase between the states \SiD 2 ) and \DiS 2 ) 
is expected to be random. This fact would in itself av- 
erage out the effect on the lifetime. In addition, if the 
superposition state is created, it is immediately (as com- 
pared to t) destroyed since, from a quantum mechanical 
point of view, we are constantly measuring the internal 
state of the ions. So any super- or subradiant effect is 
in fact expected to be destroyed, and we do not expect 
any correlation in decays from the 2 D 5 / 2 state. We have 
checked our data for correlated decays, and indeed a sta- 
tistical analysis shows no evidence of correlations. 

In Ref. mixing of the 2 D 5 / 2 state with the 2 f3/2 
state due to static electric fields was considered and found 
to be negligible. For our trap we also find this effect to 
be negligible. 

Finally, we note that the read-out time of the camera 
influences the measured duration of the shelving periods. 
As for the choice of threshold level, the read-out time has 
no effect on the measured decay rate and hence on the 
measured lifetime. 



4- The total effect of systematic errors 

Above we have identified and evaluated systematic 
errors originating from the data analysis and from de- 
shelving due to radiation and collisions. An overview of 
the estimated errors and their uncertainties is given in 
Table 13 

The effects leading to the systematic errors can be 
modelled by extra decay rates added to the natural decay 
rate: 



— - — +E 



Tnat 



or 



Tnat 



Yt 2 



,7i, 



(7) 



(8) 
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Effect 


Systematic error [ms] 


Quick re-shelving 


-1 ± 1 


Ions switching place 


-8 ±4 


De-shelving, diode lasers 


+ 1±1 


De-shelving, collisions 


+3± 1 


Total 


-5 ±4 



TABLE I: Overview of estimated systematic errors. 



where T meas is the measured lifetime as determined from 
the maximum likelihood estimate, T nat is the natural life- 
time, and the 7^'s are the extra decay rates, which can 
attain both positive and negative values in this model. 
The systematic errors given in the text and Table I cor- 
respond to T^ teas 7i. Eq. JSJ) shows that the systematic 
errors should be added linearly, yielding -5 ms, and added 
to the result of the maximum likelihood estimate, giving 
a lifetime of 1149 ms. The uncertainties of the system- 
atic errors are independent, and therefore they are added 
quadratically, yielding an uncertainty of ±4ms. 

IV. RESULT AND CONCLUSION 

By correcting the maximum likelihood estimate with 
-5 ms, as described above, we find that our final result 
for the lifetime measurement is 

Tnat = 1149 ± 14(stat.) ± 4(sys.) ms. (9) 

The largest error is the statistical, but we do have a non- 
negligible systematic uncertainty, originating from the 
correction procedure when ions change places. The asso- 
ciated error, and hence the uncertainty, could be reduced 
by increasing our signal-to-noise ratio, which would nar- 
row the time window where real events of simultaneous 
decay and shelving for two different ions could be taken 
for two ions changing place. The signal-to-noise ratio can 
be improved by frequency locking the TkSapphirc laser 
and power stabilizing the output from the doubling cav- 
ity. From Eq. we see that the statistical uncertainty 
can only be reduced by a longer data acquisition time 
(increasing N), and not simply by increasing the frame- 



rate, since the second term in the expansion is already 
negligible in our case. 

Our measurement was performed with a string of ions, 
unlike most other shelving experiments. In all single ion 
experiments, one has to consider the fact that the ion 
may heat up, so the fluorescence level drops significantly, 
and the ion can appear to be shelved, although it is not. 
In our experiment with five ions on a string, we can detect 
and discard all events of this kind since if one ion or sev- 
eral ions heat up, we would see the remaining ions move 
or heat up as well. Moreover, with a string of ions, the 
shelved ions are sympathetically cooled by the unshelved 
ions, so the number of heating events are expected to be 
reduced using a string of ions, as compared to single ion 
experiments. The major drawback of using a string of 
ions is that ions may change places, which influences the 
measured lifetime if this is not taken care of in the data 
analysis, as demonstrated here. 

In conclusion, we have measured the lifetime of the 
metastable 3d 2 D 5 / 2 state in the 40 Ca + ion using the 
shelving technique on a string of five ions. Our result 
agrees roughly at the level of one standard deviation 
with the already mentioned result obtained by Barton et 
al. 0, with the value reported in Ref. Q and with two 
theoretical papers [lfj, [lTj • On the level of two standard 
deviations, the result agrees with two other measure- 
ments where de-shelving due to diode lasers was taken 
into account 0, ^| and with the storage-ring measure- 
ment by Lidberg et al. Q • The most recent measured 
values (HQ) including our result, are in good agreement 
and group themselves around a value of about 1160 ms 
within 1% of that value. This newly obtained level of 
agreement should provide valuable input to future atomic 
structure calculations and astrononomical studies. 
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